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ABSTRACT 

We present new XMM-Newton observations of the hot-gas environments of two low- 
power twin-jet radio galaxies, 3C 66B and 3C 449, showing direct evidence for the interac- 
tions between X-ray-emitting gas and radio plasma that are thought to determine the large- 
scale radio structure of these sources. The temperatures that we measure for the two envi- 
ronments are significantly higher than those predicted by standard luminosity-temperature 
relations for clusters and groups. We show that luminosity-temperature relations for radio- 
quiet and radio-loud X-ray groups differ, in the sense that radio-source heating may operate 
in most groups containing radio galaxies. If the radio lobes are expanding subsonically, we 
find minimum ages of 3 x 10* years for 3C 66B, and 5 x 10® years for 3C 449, older than 
the values obtained from spectral ageing, which would give the radio source sufficient time to 
heat the groups to the observed temperatures for plausible values of the jet power The exter- 
nal pressures in the atmospheres of both radio galaxies are an order of magnitude higher than 
equipartition estimates of their radio-lobe pressures, confirming that the radio lobes are either 
out of equipartition or require a pressure contribution from non-radiating particles. Constraints 
from the level of X-ray emission we measure from the radio lobes allow us to conclude that 
a departure from equipartition must be in the direction of magnetic domination, and that the 
most plausible candidates for a particle contribution to lobe pressure are relativistic protons, 
an additional population of low-energy electrons, or entrained and heated thermal material. 

Key words: galaxies: active - galaxies: individual: 3C 66B - galaxies: individual: 3C 449 - 
X-rays: galaxies: clusters 



1 INTRODUCTION 

The radio lobes of twin-jet low-power (FR-I) radio galaxies ex- 
hibit a wide range of morphologies, ranging from near symmetry 
to highly irregular and distorted structure (see the 3CRR Atlas - 
|http://www.jb.man. ac.uk/atlas/ - for examples). Their large-scale 
structure is thought to be determined at least in part by interac- 
tions with the surrounding hot gas. Although earlier evidence for 
such interactions exists (Bohringer et al. 1993; Hardcastle, Wor- 
rall & Birkinshaw 1998), Chandra and XMM-Newton now provide 
the required sensitivity and spatial and spectral resolution for more 
detailed investigation, and several examples of 'holes' in intraclus- 
ter gas, with a clear correspondence to both the position and shape 
of radio-galaxy lobes, have been observed on scales of a few kpc 
with Chandra [e.g. Hydra A, McNamara et al. (2000); Perseus A, 
Fabian et al. (2000)]. Observations of the hot-gas atmospheres of 
radio galaxies have the potential to answer many questions relat- 
ing to the energetics of radio jets and lobes. One particular problem 
is how to reconcile minimum-energy radio-lobe pressures with the 
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systematically higher values determined for the external medium 
from X-ray measurements (e.g. Morganti et al. 1988; Worrall & 
Birkinshaw 2000). Study of the physical conditions and dynamics 
of the X-ray-emitting gas with good spatial resolution are important 
in clarifying this situation. 

The X-ray properties of gas in clusters are well studied (e.g. 
Sarazin 1986; Fabian 1994), and Forman, Jones & Tucker (1985) 
first showed evidence for the ubiquity of galaxy-scale hot-gas en- 
vironments in ellipticals. An analysis of the spectral and spatial 
distribution of the X-ray gas atmospheres of B2 radio galaxies was 
carried out by Worrall & Birkinshaw (2000), who found large vari- 
ation in the linear size of the atmospheres, with no correlation be- 
tween radio-galaxy size and the scale or density of the atmosphere; 
later work (Worrall & Birkinshaw 2001) found the same sample to 
have properties consistent with the luminosity-temperature relation 
for groups (Helsdon & Ponman 2000). Recent work on interactions 
between radio sources and X-ray gas has focused on the relation- 
ship between AGN and cooling-flow clusters. The details of this 
connection have become a subject of importance because of recent 
X-ray spectroscopy that conflicts with the standard cooling-flow 
models (Peterson et al. 2003; Fabian et al. 2002). Earlier work had 
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already suggested the presence of a feedback mechanism to explain 
what happens to the large quantities of cooled material in the centre 
of cooling clusters (Tabor & Binney 1993). Simulations of the evo- 
lution of radio plasma bubbles in cluster gas (Churazov et al. 2002; 
Briiggen & Kaiser 2002) show that heating can be provided in this 
way; however, the details of possible energy transfer between AGN 
and intracluster gas are still poorly understood. Although heating 
in cooling-flow regions of clusters is an important topic of study, 
the same heating processes will also occur in ordinary clusters and 
groups. Larger and more powerful FR-I sources are usually found 
in groups, where their presence is likely to have a more dramatic 
impact on the surrounding gas. 

In the more powerful FR-II radio galaxies, shock heating due 
to the supersonic expansion of radio lobes is an obvious candidate 
for a means of transferring energy from AGN to their environments. 
The first direct evidence for local heating of an X-ray atmosphere 
by radio-lobe expansion was found in Centaurus A (Kraft et al. 
2003): a rim of brighter overpressured gas capping an inner lobe has 
a temperature of 2.9 keV, while the surrounding ISM is ten times 
cooler. However, the heating process in Cen A is through shocks, as 
this inner lobe is thought to be expanding supersonically. The na- 
ture of interactions between typical FR-I radio galaxies and X-ray- 
emitting gas is somewhat different. The lobes of low-power radio 
galaxies are believed to be expanding subsonically on large scales, 
so that strong shocks would not be expected. Nevertheless, a large 
proportion of the radio galaxy's overall energy goes into displacing 
gas to inflate the lobes (Bohringer et al. 2002), and this energy must 
be dissipated into the atmosphere. This could lead to local heating 
at the lobe edges, or low-frequency sound waves may cause more 
general (and less easily detectable) heating. The presence or ab- 
sence of these effects can now be investigated with XMM-Newton. 

This paper presents the analysis of new XMM-Newton ob- 
servations of the environments of two low-power radio galaxies, 
3C 66B and 3C 449, at z = 0.0215 and z = 0.0171, respec- 
tively. These observations were obtained to look for evidence of 
interactions between the radio source and X-ray-emitting gas. Both 
sources have been extensively studied in the radio, and were cho- 
sen due to their interesting, asymmetric radio-lobe morphologies. 
3C 66B has a well-studied optical and X-ray jet (Butcher et al. 
1980; Hardcastle et al. 2001). Its radio lobes are of different shapes 
(Fig.Q right): the eastern lobe has sharp, well-defined edges, and a 
'wall' where the jet appears to stop abruptly, with some lobe emis- 
sion extending further out below the jet, whereas the western lobe 
extends into the noise. 3C 449 has similarly interesting radio mor- 
phology, shown in the detailed multi-frequency VLA observations 
of Feretti et al. (1999). It has a rounded southern lobe and narrower 
plume in the north (Fig.|2| right). i?05Ar observations (Hardcastle 
et al. 1998) showed the southern lobe to be embedded in a rim of 
hot gas, with a deficit at the position of the lobe suggesting interac- 
tions with the surrounding material, whereas the northern jet is free 
of X-ray-emitting material. 

In this paper we investigate the relationship between the X- 
ray-emitting gas and radio plasma, present new images of the ex- 
tended emission and detailed mapping of the physical conditions in 
the atmospheres, and discuss the possibility that the radio galaxies 
are heating their environments at a detectable level. 

We use a cosmology with Ho = 70 km s^^ Mpc^^, Q,m = 0.3, 
and SIa = 0.7 throughout, which gives a scale of 26.1 kpc arcmin^^ 
at the distance of 3C 66B, and a scale of 20.9 kpc arcmin"^ at the 
distance of 3C 449. 



2 DATA REDUCTION AND ANALYSIS 

We present results for the XMM-Newton EPIC MOSl, M0S2 and 
pn observations of 3C 66B and 3C 449. The duration of the 3C 
66B observation was 22227 s for MOSl and M0S2 and 19949 s 
for pn. For 3C 449 the duration was 21 109 s for MOSl, 21 1 17 s for 
M0S2, and 18498 s for pn. 

The data were reduced using the XMM-Newton Scientific 
Analysis Software (SAS) package, and filtered according to the 
methods described in the XMM-SAS Handbook. Filtering for low 
background was done by excluding time intervals with count rates 
above a threshold value obtained by studying a histogram of the 
count rate above 10 keV; the data were also filtered using the stan- 
dard flags #XMMEA_EM/P, and filtered for patterns less than or 
equal to 12 for the MOS cameras, and less than or equal to 4 for the 
pn data, as suggested in the Handbook. We additionally filtered for 
bad columns and rows. The same filtering was applied to the back- 
ground products used in our analysis (described below). Filtering 
resulted in good data of duration 17915 s, 17862 s, and 14692 s 
(MOSl, MOS2, and pn) for 3C 66B, and 20865 s, 20837 s, and 
16708 s (MOSl, M0S2, and pn) for 3C 449. 

Images and spectra were produced from the reduced data with 
the SAS task evselect, using the task evigweight to apply a vi- 
gnetting correction (as described by Amaud et al. 2002). To obtain 
the best possible images of the extended emission, several tech- 
niques were required. An interpolation routine was written that uses 
Poisson statistics to fill in the gaps due to chip boundaries and bad 
rows or columns, which would otherwise interfere with smoothing. 
The images were then adaptively smoothed using the CIAO task 
csmooth in order to facilitate identification of point sources. Con- 
taminating point sources were then removed from the original un- 
smoothed images, using the CIAO task dmfilth, which interpolates 
over point source regions using the Poisson distribution of nearby 
background regions. The resulting images were then smoothed us- 
ing Gaussian kernels to diplay the extended emission, and also 
adaptively smoothed to make more evident any compact structure 
in the extended emission. 

For each dataset, a series of spectral extraction regions were 
chosen to look for AGN components, radial and angular temper- 
ature gradients, and temperature changes related to features in the 
X-ray emission or the radio-source morphology. Figs.|3|and |4]illus- 
trate the spectral regions for 3C 66B and 3C 449 respectively. The 
extracted spectra were vignetting-corrected, as described above, 
so that on-axis response files could be used. Response files 
ml_r6_alL15.rmf, m2_r6_alL15.rmf, and epn_ff20_sdY9.rmf were 
obtained from the XMM-Newton ftp site; appropriate ancillary re- 
sponse files were generated using the SAS task arfgen. 

For spectral analysis of extended emission, it is essential to 
obtain an accurate description of the background as a function of 
position and of energy. Detailed analysis of the XMM-Newton back- 
ground by Lumb et al. (2002) has led to background template files 
useful for analysis of extended emission. Techniques have been de- 
veloped to correct for vignetted and non- vignetted components of 
the background and for variations in the soft X-ray component in 
different directions. Read (2003) has undertaken a detailed study of 
the components of the XMM-Newton background, and has created 
background event files and maps for different instrument, mode 
and filter combinations. As our observations used the medium fil- 
ter, whereas the Lumb templates (available from the XMM-Newton 
website^) are only available for the thin filter, we decided that 

^ http://xmm.vilspa.esa.es/ 
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Figure 1. Images of the 0.3-5.0 keV X-ray emission surrounding 3C 66B. The left-hand contour plot is of the point-source subtracted, smoothed (ct=22 arcsec), 
vignetting-corrected combined MOSl, MOS2 and pn data, contour levels are 0.13,0.2,0.3,0.4,0.45,0.5.0.55,0.65,0.7,0.75,0.8,0.85,0.9,0.95,1 counts/pix; the 
right hand image shows the point-source subtracted, adaptively smoothed, vignetting corrected, combined MOSl, MOS2 and pn data, with 1.4-GHz radio 
contours overlaid (from a VLA map of Hardcastle et al. 1996). Radio contour levels are \/2>2,4...512 X 1.0 X 10"'^ Jy/beam. The X-ray image scale is 2 
arcsec/pix. 




Right Ascension (J2000) 19^ Jt^m^ct 

Figure 2. Images of the 0.3-5.0 keV X-ray emission surrounding 3C 449. The left-hand contour plot is of the smoothed (cr = 32 arcsec) vignetting-corrected 
combined data; contour levels are 0.25, 0.3, 0.35, 0.4, 0.5, 0.6 , 0.7, 0.8, 1 count/pix. The right hand image shows the adaptively-smoothed, vignetting corrected, 
combined data, with 604 MHz radio contours overlaid (from a WSRT map taken from the 3CRR atlas - http://www.jb.man.ac.uk/atlas/object/3C449.htmy, 
radio contour levels are V^,2,4,..5I2 x2.0 X 10"^ Jy/beam. The X -ray image scale is 2 arcsec/pix. 
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Read's event files^ were best suited to our analysis. We filtered 
these files in the same way as our datasets, and applied the same 
vignetting correction. 

The main components of the XMM-Newton background are a 
particle component from non-X-ray events in the detector, which 
dominates at hard energies, and Galactic and extragalactic X-ray 
background, dominating at the softer end of the spectrum. The par- 
ticle background level is known to vary in intensity during and be- 
tween observations, and so in the case where a local estimate of 
background can be obtained from the source data, a scaling factor 
can be used to account for the difference in particle background 
level between source and background observations. Unfortunately 
this is not the case for either of our datasets. In both cases the ex- 
tended emission from our objects covers most of the field of view, 
so that it was not possible to obtain an estimate of the local back- 
ground level. This is particularly problematic in the case of 3C 449, 
where a large molecular cloud in our Galaxy lies along the line-of- 
sight, so that the absorbing column density is high. In both cases 
scaling the background spectra would potentially introduce more 
uncertainty into an already poorly determined component of the 
background. 

We determined the necessary scaling factors between the par- 
ticle levels in our observations and the background files by using 
counts measured in the out of field-of-view regions (due entirely to 
non-X-ray events) in our energy range of interest, 0.3 - 5.0 keV, 
as we found that there was significant spectral variation in particle 
background level, but little spatial variation across the detector (be- 
fore applying the vignetting correction). To determine the amount 
of error introduced into our spectral analysis by our background, 
we used these scaling factors to determine the expected contribu- 
tion of particle background to the count level in the region of low- 
est surface brightness for each dataset (an outer annulus), where the 
background contribution would be most important. We found max- 
imum offsets in 0.3-5.0 keV counts of 0.5, 3 and 7 per cent (MOS 1, 
M0S2, and pn) for 3C 66B, and 9, 8 and 16 per cent for 3C 449, 
between images with scaled and unsealed background subtraction. 
As this uncertainty in counts is not negligible, particularly in the 
case of 3C 449, we compared the mekal model fits obtained using 
scaled and unsealed background spectra, and found that the best- 
fitting temperatures did not differ significantly for the worst case 
scenarios for either 3C 66B or 3C 449. There was a small system- 
atic difference in the best-fitting abundances for the pn spectrum of 
3C 449; however, the unsealed value is more physically plausible 
and in better agreement with the other two cameras. We concluded, 
therefore, that the effect of background scaling on the measured 
physical parameters is largely insignificant, particularly compared 
to the unknown error due to the difference in soft Galactic back- 
ground between our source and background datasets, and so we 
used unsealed background data in our XSPEC spectral fitting. 

In order to study the spatial distribution of the X-ray-emitting 
gas, we obtained radial surface-brightness profiles of the two ob- 
jects. The profiles were fitted with convolved point-source- and 13- 
models following the method of Birkinshaw and Worrall (1993), 
individually for each camera, so as to correctly model the differ- 
ent camera responses. The PSF components were modelled using 
the analytical description provided in the document XMM-SOC- 
CAL-TN-0022 obtained from the XMM website. The background 
determination for this stage of the analysis was done in a different 
way to that of the spectral analysis, as a small error in the cho- 

^ available from |http://www.sr.bham.ac.uk/xmm3/BGproducts.htmI| 



sen background level could have a significant effect on the shape 
of the radial profile; in particular, the contribution of the particle 
background is a function of radius, since the vignetting correction 
is incorrectly applied to the particle-background component in our 
approach. Therefore, for this part of the analysis, we used a double- 
subtraction technique. We determined background levels in each 
annulus from the Read background files; these were then scaled by 
the factors given above and subtracted from the source counts in 
each bin. We then subtracted the remaining count density in an out- 
ermost annulus (already background corrected in the first stage of 
background subtraction) from all the bins, so as to correct for any 
differences in diffuse Galactic background level between the source 
and background datasets. 

Once temperatures had been determined from spectral fitting, 
and the surface-brightness profiles modelled as described above, 
this information was used to determine densities and pressures in 
the X-ray-emitting gas following the method of Birkinshaw and 
Worrall (1993), to compare these with the properties of the radio 
lobes. 

3 RESULTS 
3.1 3C 66B 

Fig.Qshows a contour plot of a smoothed image obtained by com- 
bining data from the three XMM-Newton cameras, as well as an 
adaptively smoothed combined image with radio contours overlaid. 
The morphology of the diffuse emission surrounding 3C 66B has 
a number of interesting features: there appear to be two deficits 
in surface brightness to the west and east of the centre, and there 
is a region of brighter emission at the edge of the eastern deficit. 
We tested the significance of the deficits by comparing the counts 
in these regions and in non-deficit regions at the same radius, and 
found the eastern deficit to be significant at the 4cr level, and the 
western deficit to be significant at the la level. The overlaid radio 
contours on the adaptively smoothed image illustrate the relative 
positions of X-ray and radio features. There is a correspondence 
between the positions of the X-ray deficits and the radio lobes, and, 
in addition, there is a brighter blob of emission at the end of the 
eastern jet/lobe. The relationship between these radio and X-ray 
features is discussed in more detail in Section 4.1 below. 

3C 66B has been observed in the X-ray with Chandra (Hard- 
castle, Birkinshaw Sl Worrall 2001) and the observed X-ray emis- 
sion was found to comprise an AGN component, emission from 
a jet, and extended emission attributed to hot gas extending to 
at least 40 arcsec from the core. The resolution of XMM-Newton 
means that the AGN and jet are not spatially separable from the 
extended emission in our data, and so we obtained a spectrum for 
a region consisting of the central 60 arcsec, to include all the AGN 
emission, for comparison with the Chandra results. A background 
spectrum was taken from a surrounding annulus to facilitate com- 
parison (in contrast to the background method described above and 
used for all other spectra presented here). We fitted the data from 
the three XMM-Newton cameras with an absorbed power law plus 
mekal model, and found a best-fitting power-law energy index of 
O-S^laiT' consistent with the value of 1.14±0.05 of Hardcastle et 
al. (2001) , and best-fitting temperature for the extended component 
of 0.65±0.05 keV, also in agreement with their results (all errors 
quoted for our analysis are Icr confidence levels; those for mekal 
temperatures and abundances are Icr for two interesting parame- 
ters). The model used the Galactic Nh value of 8.36 x 10^" cm~^ 
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Figure 3. Spectral extraction regions for 3C 66B; labels indicate the regions listed in TablelTI except dl and d2 which are the regions used for analysis of the 
nearby cluster and 3C 66A respectively. Region a4 excludes dl and d2 from the extraction region, and regions cl,2 and 3 exclude a circle of radius 60 arcsec 
centred on the nucleus. 



(Murphy et al. 1996), gave a best-fitting abundaiice of O-OSlJi 
solar, and had a reduced of 0.98 (197 d.o.f.). We extracted a 
matching spectrum from the same region of the Chandra data, and 
found that a joint fit to all four spectra (MOSl, M0S2, pn, and 
Chandra) resulted in fits which did not differ significantly from 
the results quoted above. The high degree of consistency between 
the XMM-Newton and Chandra spectra suggests that there has not 
been significant variation of the AGN between the observations. 
The Chandra data better constrain the AGN spectral parameters 
because a single power-law model could be applied to the spa- 
tially separated spectrum. We measure a total unabsorbed flux (in 
the energy range 0.5 to 4.0 keV) from the power-law component 
(1.90±0.47) X 10"" ergs cm"^ s"\ consistent with the total flux 
of core and jet components from Hardcastle et al.'s Chandra results 
of 2.31 xlO"" ergscm"^ s"^ 

A spectrum was obtained for the entire extent of the diffuse 
X-ray emission surrounding 3C 66B out to a radius of 450 arcsec, 
excluding the central 60 arcsec. The best-fitting model has a tem- 
perature of 1.73lo:°4 keV with an abundance of 0.28±0.03 ix^ /n 
= 1.14 for 718 d.o.f.). The total X-ray luminosity (in the range 0.5 
to 10 keV) was found to be (2.81±0.05) x 10^^ W. We applied 
the X-ray luminosity-temperature relation for groups described in 



Helsdon & Ponman (2000) [found by Worrall & Birkinshaw (2001) 
to be consistent with a sample of FR-I radio-galaxy environments], 
which gives an expected temperature of 0.79±0.03 keV, signifi- 
cantly lower than that obtained from our observations. We return to 
this point in Section 4.3. 

As described in Section 2 above, a series of spectral extraction 
regions were chosen for more detailed analysis; these included an- 
nuli for radial temperature analysis, pie-slice regions for study of 
angular variations, and regions selected to study features of partic- 
ular interest evident in the image of Fig.Q The choice of extraction 
regions is illustrated in Fig.|3| The regions used for study of two 
other objects in the field of view, 3C 66A and another unidentified 
extended object (discussed in Sections 3.1.1 and 3.1.2), are also 
shown in this Figure. The results of spectral analysis are given in 
TableQ which shows the results of jointly fitting an absorbed mekal 
model to spectra from the MOSl, M0S2 and pn cameras in the en- 
ergy range of 0.3 to 5.0 keV. The neutral hydrogen column density 
was fixed at 8.36 x lO*^" cm"^ (Murphy et al. 1996). Temperature 
and abundance were allowed to vary. Best-fitting models are shown 
in Table Q and plots of temperature as a function of angle and ra- 
dius are shown in Fig.|5| Fig.Qshows the best-fitting mekal model 
for the total emission. 
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Figure 4. Spectral extraction regions for 3C 449: labels indicate the regions listed in Table I2I Region cl excludes a circle of radius 60 arcsec centred on the 
nucleus. 



Radial surface brightness profiles were extracted for each 
camera and fitted with convolved point-source plus /3-models as 
described in Section 2 above. The values for the three cameras 
were combined for each fit to determine the best joint fit to all the 
three profiles. We found a best-fitting model with /3 = 0.375lo , 
core radius = 179.6t^^ arcsec with x^/n of 1.88 (267 d.o.f.) (er- 
rors are la for two interesting parameters). An unacceptable fit is 
not surprising, as the images and spectral analysis show that there is 
clearly small-scale structure in the extended emission. Fig.|6|shows 
the pn radial profile with the best-fitting model to the individual 
profiles. There are systematically positive residuals in the fits at 
small radii, which could be due to a galaxy-scale atmosphere. Our 
spectral results, and the Chandra data of Hardcastle et al. (2001) 



indicate that a galaxy-scale atmosphere would be expected in ad- 
dition to the larger group-scale gas fitted by the model described 
above. We fitted a second /3-model component by holding the large 
/3-model fixed, and found that including a model with (3 = 0.9 
and core radius of 20 arcsec significantly improved the fit for both 
MOSl and MOS2 (though the fit is still unacceptable for the rea- 
sons explained above), and resulted in a slight, but not significant 
improvement for the pn fit. The two-/3-model fit for MOS 1 is shown 
in Fig.|6| The presence of this smaller-scale component is consis- 
tent with the Chandra data; however, its parameters cannot be well- 
constrained by our XMM-Newton data, where much of the emission 
on these scales is hidden by the PSF, or by the Chandra data, since 
its surface brightness is rather low. 



XMM-Newton observations of the hot-gas atmospheres of 3C 66B and 3C 449 7 



Region Net counts Temp (keV) Abundance /n (d.o.f.) 



tot 

al 

a2 

a3 

a4 

bl 

b2 

b3 

b4 

b5 

b6 

b7 

b8 

b9 

blO 

bll 

bl2 

cl 

c2 

c3 

c4 



37593 

7029 

9070 

8748 

12613 

1235 

1392 

1359 

1235 

1138 

1478 

1250 

1385 

1056 

1562 

1181 

1101 

900 

2260 

5120 

1401 



3+0.03 
-0.04 

1.99±0.11 

f-O.lO 
-0.12 
^+0.09 
-0.06 
^+0.05 
-0.02 
7+0.28 
-0.19 

q + 0.45 

-0.25 
7+0.24 
-0.25 
.;+0.56 
-0.43 
D+0.49 
'-0.33 
,+0.43 
-0.31 
3+0.33 
-0.25 
2 30+"-''° 
1+0.32 
-0.21 
-1+0.31 
-0.28 
^+0.32 
-0.25 
f-0.22 
-0.09 
2+0.34 
^-0.26 
9 ,9+0.20 
^■^^ -0.22 



1.73^ 



2.043 

1. soj 

1.44j 
1.97^ 
2.25^ 
2.07^ 
2.55^ 

2. ( 
2.23^ 
1.99^ 



1.81] 
2.101 
1.96j 
1.48j 
1.; 



2.05 
2.36 



+0.14 
0.15 
+0.46 
0.34 



0.28±0.03 
0.45j 
0.40^ 



+0.09 
0.08 
+0.08 
0.07 

0.31 ±0.06 

+0.03 
0.02 

0.64^ 



0.18 

1+0:39 

-0.24 
7+0.47 
-0.23 
1+0.35 
-0.22 
-1+0.43 
-0.25 
44+0.28 
"■^^-0.19 
7+0.26 
-0.17 
3+0.22 
-0.13 
,+0.35 
-0.25 
7+0.19 
-0.13 
S+0.21 
-0.15 
S+0.28 
-0.16 
,+0.14 
-0.09 

0.42 



0.67: 
0.64j 
0.50^ 



0.37: 
0.291 
0.631 
0.271 
0.381 
0.381 
0.23^ 



-,+0.25 
-0.20 



1.14(718) 
1.02 (269) 
1.15 (330) 

1.17 (311) 
1.02 (395) 
0.90 (54) 
0.97 (61) 
0.93 (59) 
1.21 (54) 
1.70 (50) 
1.37 (65) 
1.13 (56) 
1.25 (62) 
1.00(46) 
0.94 (69) 
1.20(51) 
0.89 (48) 
1.27 (39) 
1.43 (97) 
0.96 (204) 

1.18 (62) 



Table 1. 3C 66B: Best-fitting mekal models; the choice of regions is described in Section 3.1 and illustrated in Fig. |3] Errors are Icr for two interesting 
parameters. 
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Figure 5. Plots of temperature variations in 3C 66B's environment as a function of radius and angle with abundace for the fits shown in Table l3.l1 



3.1.1 3C66A 

3C 66A is a BL Lac object at a redshift of 0.44, seen 6.5 arcmin to 
the north-west of 3C 66B in the images of Fig.Q Initially thought 
to be related to 3C 66B, it was soon determined to be a separate 
object (Mackay et al. 1971; Northover 1973), and was optically 
identified by Wills & Wills (1974). Its spectrum was extracted for 
the three cameras, and an absorbed power law model was fitted 
in the energy range of 0.3 to 7.0 keV. Allowing the Galactic Nh 
to vary, a photon index of 2.47 ± 0.03 was obtained for a value 
of Nh = (6.9±0.5) x 10^" cm"^ (x^/n = 1.17 for 668 d.o.f.). 
This value for Nh is slightly lower than the value of 8.36 x 10^" 
cm~^ used in the 3C 66B spectral fits above (Murphy et al. 1996). 



Fixing the Nh at the Galactic value, a photon index of 2.55±0.01 
was obtained (x^ /n = L19 for 669 d.o.f.). The photon indices for 
both fits are consistent with the results of Sambruna et al. (1994), 
Fossati et al. (1998), and Worrall & Wilkes (1990) [as tabulated by 
Donato et al. (2001)]. The main residual component in these fits 
was at high energies; including a flat or heavily absorbed second 
power-law component significantly improves the fit. For a second 
power law with photon index of -0.3±0.5, we obtain x^/*^ of 1.1 1 
(667 d.o.f.); a heavily absorbed power law (Nh = 22^^* x 10^^ 
cm^^) with energy index 2.11^ 3 gives x^/^ - 1-10 (666 d.o.f). 
The photon index and normalization of the single component fit 
with fixed Galactic column density were adopted, and this model 
led to a 1-keV flux density of 0.378 /^Jy, which is lower than the 
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Figure 6. a) Radial surface brightness profile for 3C 66B using pn data only; model is the best-fitting point-source plus /3-model to the MOSl, M0S2 and pn 
profiles, with /3 = 0.375 and Xc = 179.6 arcsec, b) MOSl profile with the same model as for a), but with a second /3-model component having {3 = 0.9 and tc = 
20 arcsec 




values obtained in the three earlier studies by a factor of 3 to 5. 
This amount of variabihty is not implausible, as optical monitoring 
shows variability of 1.5-2 mag (e.g. Takalo et al. 1996). There was 
no significant variability evident in the light-curve of 3C 66A over 
the course of our observation. We fitted a radial surface-brightness 
profile to the emission from 3C 66A, and this was adequately fit by 
a point-source model with no extended emission. Although this is 
a bright object, the count rates are too low for pile-up to occur, as 
is confirmed by SAS task epatplot, so we conclude that the lower 
X-ray flux density is a real effect of source variability. 



3.1.2 Nearby cluster 

The images of Fig. Q show an extended object ~ 7 arcmin to the 
south-east of 3C 66B. We find 4874 background-subtracted counts 
(MOSl, M0S2 -I- pn) in a 78 arcsec circle. The position of this 



emission does not coincide with any previously detected cluster. 
An absorbed mekal model was fitted in the energy range 0.3 to 7.0 
keV, with an abundance of 0.4 solar, allowing redshift and tempera- 
ture to vary. The best fit temperature was 3.7lo3keV, ataredshift 
of 0.35l[;;J!4; xV"- for this fit was 1.33 for 174 d.o.f. The unab- 
sorbed X-ray luminosity (between 2 and 10 keV) was found to be 
(1.66^009) X 10''' W. Applying the X-ray luminosity-temperature 
relation (Amaud & Evrard 1999) results in an expected tempera- 
ture of 3.0 keV which is in approximate agreement with our best 
fit value. The X-ray luminosity/brightest cluster member relation 
of Edge & Stewart (1991) predicts that the brightest cluster galaxy 
should have a visual magnitude of ~ 19. DSS images of this field 
show a faint extended object at the position of the centre of emis- 
sion, with several other objects of similar magnitude nearby. There- 
fore we conclude that this is a strong candidate for a new X-ray- 
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detected cluster of galaxies at z — 0.35 at a = 02:23:34, S = 
+42:54. 



3.2 3C 449 

Fig.|5|shows a contour plot of a smoothed image obtained by com- 
bining data from the three XMM-Newton cameras, as well as an 
adaptively smoothed combined image with radio contours overlaid. 
The radio-X-ray overlay illustrates the presence of a deficit at the 
position of the southern radio lobe; we tested the significance of 
this deficit by comparing the count rate with that of a region at sim- 
ilar radius in the direction perpendicular to the jet axis, and found 
it to be significant at greater than 5a. 

Hardcastle et al. (1998) presented RO SAT imaging and spec- 
tral analysis of the gas surrounding 3C 449. We took a spectrum of 
the full extent of the X-ray emission out to a radius of 420 arcsec 
to compare with the ROSAT results. We found a best-fitting mekal 
model with temperature of 0.98±0.02 keV, consistent with their 
results at the 2a level, and a luminosity in the energy range 0.2 
- 1.9 keV of (1.73±0.05) x 10^^ W, consistent with their value 
of 2.1 X 10'^'' W (corrected for the different cosmology we use 
here). In the energy range 0.5 to 10 keV, we measure a luminos- 
ity of (1.98±0.06) X 10^^ W. The luminosity-temperature relation 
of Helsdon & Ponman (2000) predicts a temperature of 0.73±0.05 
keV, lower than the value we measure. We took a spectrum of the 
central 60 arcsec to look for an AGN component; however, the best- 
fitting mekal plus power-law model has a photon index of 0.63lQ g2 
(with a temperature of 1.37i;5:w and /n of 1.26 for 93 d.o.f). 
If the photon index is held fixed at the more physically plausible 
value of 2 (i.e. an energy index of I), a x^/n of 1.27 (for 94 d.o.f.) 
is obtained. This is not a significantly worse fit; therefore we can- 
not constrain usefully the power-law index with these data. Taking 
the second model with fixed power-law index, we find a flux from 
AGN and jet components of 7.7 xlO^^* ergs cm~^ s^^, a factor 
of 3 lower than that of 3C 66B. This is consistent with the flux in 
the point-source component of our radial profiles described below. 

We defined a set of similar spectral extraction regions to those 
used for 3C 66B to look for radial and angular temperature struc- 
ture. The choice of regions is illustrated in Fig.|4| Fits were per- 
formed in the same way as for 3C 66B, using a fixed Galactic col- 
umn density of 11.8 x 10^" cm"^ (Murphy et al. 1996): we found 
that varying the column density did not significantly improve our 
fits. Table |2| shows the best-fitting mekal model fits for fixed and 
free abundance. Plots of temperature variations, shown in Fig.[S|il- 
lustrate that the temperature is lower at position angles of around 
180 and 360 degrees, corresponding to the position of the radio 
jets and lobes, whereas higher temperatures are found where there 
is brighter X-ray gas, perpendicular to the radio axis. The tempera- 
ture is also found to decrease with radius, from 1.39 keV at a radius 
of 45 arcsec to 1.04 keV at a radius of 260 arcsec. 

Radial surface brightness profiles were extracted for each 
camera and fitted with convolved point-source plus /3-models as for 
3C 66B. We found a joint best-fitting model with /3 = 0.375±0.025, 
core radius = 48.4l^o''5 arcsec with /n of 1.27 (267 d.o.f.). As 
in 3C 66B, the unacceptable fit is likely to be due to the dumpi- 
ness of the environment. Fig.|9|shows the pn radial profile with the 
combined best-fitting model. The radial profiles show no evidence 
for a galaxy-scale atmosphere. 




300 



Figure 9. Radial surface brightness profile for 3C 449 using pn data only; 
the model is the best-fitting point-source plus /3-model to the MOS 1 , MOS2 
and pn profiles, with /3 = 0.375 and rc = 48.4 arcsec 



4 DISCUSSION 

4.1 Interactions of radio plasma and X-ray-emitting gas 

4.1.1 3C66B 

The contour plot and adaptively smoothed image of3C66B (Fig.[n 
show strong evidence for interactions between a low-power (FR-I) 
radio galaxy and its environment. The presence of deficits in the 
X-ray surface brightness coincident with the extended radio struc- 
ture provides evidence both for the influence of the X-ray gas on 
the lobes, and for the displacement of X-ray-emitting gas. Hot gas 
appears to surround the eastern radio lobe, which has a sharp edge 
and 'wall' at the end of the jet; in contrast, the more diffuse western 
lobe is relatively free of gas. This suggests that the eastern lobe has 
been constrained into its present shape by the surrounding hot gas, 
whereas the absence of gas at similar distances from the core on 
the opposite side of the source left the western jet relatively free to 
propagate. That the eastern lobe is the shorter of the two by a factor 
of ~2.5 is consistent with the suggestion that the environment has 
slowed the lobe expansion. Emission round the edges of the west- 
em lobe for a short distance indicates that the expanding lobe has 
pushed material out of the way; however, the lobe morphology sug- 
gests that the gas was originally less dense than on the eastern side, 
allowing the western lobe to expand faster than its counterpart. 

A particularly interesting feature in 3C 66B's X-ray-emitting 
environment is the bright blob of gas at the end of the eastern jet, 
just beyond the sharp boundary in the radio emission. As shown 
in Fig.Q the radio jet appears to end suddenly, and lobe emission 
appears below and to the south of the blob of gas. However, it is 
unclear whether the blob of gas has been pushed out from near the 
centre of the radio galaxy, or whether the lobe has only recently 
encountered a denser clump of gas, which we find to have a mass 
of 4 X lO'' Mq (obtained from the volume emission measure of 
the best-fitting spectrum). 

Without detailed hydrodynamical simulations it is difficult to 
determine the dynamics of this gas blob in the atmosphere. If it 
originated close to the centre of 3C 66B, we can determine a lower 
limit to the time required for the jet to push the blob out to its cur- 
rent distance from near the centre by neglecting the effect of grav- 
ity, and assuming that all of the underlying jet's thrust acts on the 
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Table 2. 3C 449: Best-fitting mekal models; the choice of regions is described in Section 3.2 and illustrated in Fig.l4l 
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Figure 8. Plots of temperature variations in 3C 449's environment as a function of radius and angle for the fits shown in Tablel2l 



blob. Assuming a jet power of 3.5 x 10 W [scaling the jet power 
of 3C 31 from models of Laing & Bridle (2002) by the ratio of 
low-frequency radio luminosity in the two objects], we find that it 
would take 1.4 x 10* years. This is likely to be less than the age of 
the source, as discussed in Section 4.3. However, gravity is likely 
to significantly increase the time required. 

We calculated the mass swept out by the front of the expand- 
ing radio lobe to be ~4 x 10^ ^0' the same as the mass we calcu- 
late for the blob. This mass was determined by integrating the den- 
sity of the /3-model atmosphere, assuming the volume swept out by 
the lobe is a cone of half-opening angle 45°. We also find that the 
temperature increase produced by adiabatic compression of the gas 
evacuated from the lobe to the current blob volume is consistent at 
the 2a level with our measured temperature of 2.36lg'34 keV (sig- 
nificantly hotter than the surrounding medium) assuming 7 = 5/3. 
However, it seems unlikely that this could explain its origin, as the 



gas would be expected to spread out, rather than form such a dense 
compact blob. 

The position of this blob and the radio-source morphology at 
the end of the lobe are suggestive of a scenario where the gas forms 
an obstacle choking off the flow and forcing the lobe to expand 
laterally. Either the jet could have encountered a denser clump of 
gas which is part of the environment, or an infalling clump of gas. 
There is no evidence that the blob has a different redshift than the 
surrounding material. Both scenarios are physically plausible and 
consistent with the radio-source morphology of 3C 66B. 

4.1.2 3C449 

Evidence was presented in Hardcastle et al. (1998) for a deficit in 
X-ray surface brightness coincident with the southern radio lobe of 
3C 449 and Fig. |2| confirms the presence of this deficit. The nar- 
row northern lobe is free of X-ray-emitting gas at 100 kpc from the 
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core, whereas there is more gas to the sides and below the more 
rounded southern lobe (unfortunately the southern lobe extends to 
the edge of the field of view, so the amount of gas below the radio 
lobe is not clear). In both cases there are deficits in X-ray surface 
brightness at the lobe positions. The 1.4-GHz radio map of Fer- 
etti et al. (1999) shows the southern lobe to have more internal 
structure than appears in our lower resolution map. However, the 
rounded appearance of this lobe suggests a similar scenario to that 
described above for 3C 66B: the southern lobe expansion along the 
jet axis has been slowed by the denser material we see producing 
the higher level of X-ray emission in that direction. Therefore, our 
images of the X-ray-emitting gas surrounding 3C 449 give further 
support for the idea that the X-ray environment plays an important 
role in determining the morphology of the radio lobes of low-power 
radio galaxies. The radio maps of Feretti et al. (1999) also show 
smaller "bulges" in both jets nearer to the core, which could be due 
to the jets crossing the boundary between a denser galaxy atmo- 
sphere and the group gas. We find that the gas in the surrounding 
rim of 3C 449's southern lobe has a mass of 1.1 xlO^° Mq (we 
took the spectrum from a section of the rim using local background 
subtraction, then extrapolated this mass to a sphere of thickness 34 
arcsec; this is likely to be an overestimate, since there will be no 
rim on the core side of the lobe). As shown in Table|2| the rim has 
a temperature of 1.01 ±0.03 keV, consistent with the temperature in 
the outer regions of the extended gas. We calculate that the amount 
of gas swept up by the lobe front is ~ 6 x 10^" Mq. It is therefore 
plausible that most of the rim gas has been evacuated from the lobe 
and moved to its current position by lobe expansion. 

4.2 Pressure balance between radio lobes and X-ray 
environment 

The radio lobes of low-power radio galaxies are thought to be ex- 
panding subsonically on the largest scales, so that the internal pres- 
sure within the radio lobes should be similar to, but slightly higher 
than that of the external environment. However, internal radio-lobe 
pressures obtained by assuming an electron filling factor of unity, 
equipartition of energy in particles and magnetic fields, and that 
the only contribution to pressure comes from the population of 
synchrotron-emitting electrons, are found to be significantly lower 
than the external pressures obtained from X-ray measurements (e.g. 
Morganti et al. 1988; Worrall et al. 1995; Worrall & Birkinshaw 
2000). If the lobes are expanding supersonically, these problems 
are exacerbated (more detailed discussion of expansion speeds fol- 
lows in Section 4.3). We determined radial pressure profiles in the 
X-ray enviroiraients of 3C 66B and 3C 449; these are shown in 
Fig. II II Internal radio-lobe pressures were determined using the 
code of Hardcastle et al. (1998), choosing an electron energy spec- 
trum with a power-law number index of 2, minimum energy of 5 
X 10® eV and maximum energy of 6 x 10^^ eV. 

For 3C 66B, we find the external pressure acting on the eastern 
lobe of 3C 66B at a radius of 75 arcsec to be a factor of ~24 higher 
than the equipartition internal lobe pressure, and for the southern 
lobe of 3C 449 we find the external pressure at a distance of 350 
arcsec to be ~16 above the internal lobe pressure. The result for 
3C 449 agrees with Hardcastle et al.'s ROSAT results, and both re- 
sults are consistent with ROSAT measurements for similar sources 
(e.g. Worrall et al. 1995). Therefore, as the radio lobes cannot be 
underpressured, we must consider which of the assumptions used 
to determine the internal lobe pressures are incorrect. 

In order for the population of synchrotron-emitting electrons 
to provide the necessary pressure without a significant contribution 
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Figure 10. Plot of the modelled spectium of 3C 449's southern lobe for 
magnetic field strength of 20 pT (to achieve pressure balance with external 
environment). The solid line is the synchrotron emission, passing through a 
radio data point. The second data point is an upper limit to X-ray inverse- 
Compton emission as described in the text, which is significantly below the 
dashed line indicating the expected inverse Compton emission from CMB 
photons scattered by the population of synchrotron-emitting eletrons. The 
lower dashed line is the expected level of synchrotron self-Compton emis- 
sion 



from other particles, conditions in the radio lobes of the two ob- 
jects must deviate substantially from equipartition. For 3C 66B, we 
find an equipartition magnetic field strength of 0.34 nT (1 nT = 
lO/iG). To obtain pressure balance, field strengths of 30 pT or 2.0 
nT are required. In 3C 449 we obtain similar results: the equipar- 
tition magnetic field is 0.23 nT, whereas field strengths of 20 pT 
or 1.2 nT are required for pressure balance. However, from spec- 
tral fitting to the X-ray deficit regions coincident with these radio 
lobes we obtain an upper limit on the flux from inverse Compton 
emission of 10 nJy (3C 66B) and 3 nJy (3C 449), which leads to 
lower limits on magnetic field strength of 90 pT (3C 66B) and 80 
pT (3C 449). With such fields, the lobes would still be underpres- 
sured by factors of 5 (3C 66B) and 7 (3C 449). Fig.[lOlshows that 
the level of inverse Compton emission expected from the lobes of 
3C 449 for a field strength of 20 pT (required for pressure balance 
in the particle-dominated regime) is significantly higher than the 
upper limit of a power-law component of X-ray flux, as taken from 
fitting a mekal plus power-law model to the spectrum from the lobe 
region. Therefore we conclude that in the scenario where no other 
particles are present, if the lobes of these two objects deviate from 
equipartition conditions, it must be in the direction of magnetic 
domination, requiring field strengths of ~ 1-2 nT. Such magnetic 
domination is not found where direct measurements of magnetic 
field strengths have been made from X-ray inverse Compton emis- 
sion. In no case is the field strength significantly higher than the 
equipartition value; Feigelson et al. (1995) found a level of X-ray 
inverse Compton in the lobes of Fornax A slightly higher than, but 
consistent with equipartition, and Tashiro et al. (1995) found the 
lobes of PKS 1343-601 to be particle dominated. 

Our results also constrain pressure contributions from non- 
radiating particles. Although the presence of some entrained ma- 
terial in the lobes is necessary (e.g. from models of jet decelera- 
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tion), the X-ray deficits rule out a significant contribution to pres- 
sure from a thermal gas component at the temperature of the sur- 
rounding atmosphere. The presence of sufficient cold gas would be 
expected to produce significant Faraday depolarization of the lobes, 
which is not observed (e.g. Feretti et al. 1999), although it is pos- 
sible to devise magnetic field structures where such material could 
be present but does not result in significant depolarization (Laing 
1984). However, it is implausible that large quantities of material 
could be hidden in this way. Entrained material could be heated 
(by the entrainment process, or by energy transfer from interactions 
with the relativistic electrons), and so could provide sufficient pres- 
sure without producing detectable X-ray emission: e.g. in 3C 66B, 
the missing pressure might be provided by gas at -^5 keV. Missing 
pressure could also be hidden in an extra population of low-energy 
electrons, producing inverse Compton emission at a level we could 
not detect. However, these electrons must have 7 ~ 1 — 200, requir- 
ing the power-law number index to steepen to 3 below ~ 10* eV, 
then to drop off sharply at energies of ~ 5 x 10^ eV, since large pop- 
ulations of electrons at these energies would also be expected to re- 
sult in depolarization. This model is somewhat implausible. There 
could also be a contribution from relativistic protons, although this 
would require large proton to electron ratios of k ~260 for 3C 
66B, and ~130 for 3C 449. It also possible that there is a plasma 
filling factor 4> less than unity; however, to achieve pressure bal- 
ance, <7i - 1/260 (3C 66B) and 1/130 (3C 449) would be required. 
More detailed discussion of this scenario is found in Hardcastle et 
al. (2000). 

We have assumed that the source is in the plane of the sky, 
which may not be the case. As described in Hardcastle et al. (1998), 
if the radio source is at an angle of 6 to the line of sight, the mini- 
mum internal pressure falls as (sin 9)'^^^ , so that small angles to the 
line of sight are required to achieve pressure balance. In 3C 66B, an 
angle of 7° to the line of sight is required, and in 3C 449 an angle of 
11° is required. In both cases these angles are implausible because 
they imply extremely large physical sizes. Additionally, Feretti et 
al. (1999) find that in 3C 449, 6 > 75°, which argues strongly 
against such an explanation. 

4.3 Radio-source heating of the X-ray environment 

In Section 3.1 we showed that the environments of 3C 449 and 
particularly 3C 66B are significantly hotter than predicted temper- 
atures from the luminosity-temperature relation of Helsdon & Pon- 
man (2000). It is interesting to speculate on possible explanations 
for these increased temperatures. Expanding radio lobes would be 
expected to put a significant proportion of their energy budget into 
PAV work on their environment. Only one example of radio-source 
heating of the environment has been directly detected so far: shock- 
heating of the environment of Cen A (Kraft et al. 2003). In the 
case of 3C 66B, determining the atmosphere's heat capacity using 
AE/AT = (3/2)nfc, where n = 2.18np, which is the total num- 
ber of particles, where rip is the number of protons (determined 
from the volume emission measure obtained from spectral fitting 
to the entire environment to a radius of 450 arcsec), we calculate 
that the energy needed to heat the entire environment by 0.9 keV 
(from the predicted temperature of ~0.8 keV to the observed 1.7 
keV) is ~ 2.4 x 10^^ J. If we assume that the radio jets of 3C 66B 
have a luminosity of 3.5 x lO^'^ W [scaling the jet luminosity of 
3C 31 (Laing & Bridle 2002) by the ratio of low-frequency radio 
luminosities of the two objects] and that 1/3 of this energy goes 
into heating the environment by PdV work (since pressure is 1/3 
X energy density), we find that it would take ~ 3 x 10** years for 



3C 66B to heat its environment by this amount. Using 1.4- and 8- 
GHz radio data, we fitted a Jaffe & Perola (1973) aged synchrotron 
spectrum with injection index of 0.5, in a magnetic field of 0.3 nT, 
to material at the far edge of the eastern lobe, we obtain a spectral 
age for 3C 66B of ~ 10** years. In order to provide the necessary 
energy input for a temperature increase of 0.9 keV in 10* years, the 
jets would need to be ~ 3 times more powerful than the value we 
use above (or a larger proportion of the jet energy would have to go 
into radio-source heating). 

However, spectral ages are based on the spectrum of the oldest 
measured radio material in the lobes. Such ages should be consid- 
ered as lower limits, as there may be even older material which 
is no longer emitting. To investigate whether 3C 66B and 3C 449 
could be older than their measured spectral ages, we used the sizes 
of the two radio sources to estimate their ages, assuming they are 
in the plane of the sky (thus providing a lower limit to size and 
hence age). The distance from the radio core to the western lobe 
edge of 3C 66B (in Fig.0 is 175 kpc, and the sound speed in its 
environment is ~ 7 x 10^ m s~^, which means that if the lobe 
front is moving subsonically, the source age must be at least 3 
xlO* years. In the case of 3C 449, the northern edge (in Fig.|2j 
is 290 kpc from the core and the sound speed is ~ 6 x 10° m s~^, 
which implies a source age exceeding 5 x 10* years. Interestingly, 
in both cases these values are close to the times required to raise the 
group temperatures from the predicted radio-quiet temperatures to 
those observed. There are two reasons to favour radio-source heat- 
ing as an explanation of hotter group gas: firstly, in both sources the 
morphology suggests that the lobes are unlikely to be supersonic; 
secondly, were the lobes supersonic, they would be expected to be 
overpressured, increasing the need for additional pressure compo- 
nents or magnetic field as described in Section 4.2. Therefore we 
suggest that FR-I sources are likely to be older than typical spectral 
ages suggested by earlier work (e.g. Andernach et al. 1992; Parma 
et al. 2002). This means that both 3C 66B and 3C 449 could be 
old enough to have heated their environments by the amount nec- 
essary to produce the observed temperature increases. Eilek (1996) 
reached similar conclusions about radio-source ages from a more 
detailed dynamical analysis. 

The luminosity-temperature relationship we have used is 
based on an extrapolation of the relation for clusters to radio-galaxy 
environments (Worrall & Birkinshaw 2000). If radio-source heating 
is occurring in the environment of 3C 66B, we might expect that it 
could affect the temperatures of other environments in their sample, 
so that we are not comparing 3C 66B with unheated environments. 
The larger X-ray group sample of Helsdon & Ponman (2000) also 
contains a significant proportion of radio-galaxy host groups. They 
find a relation for groups which differs from the cluster relation. 
Therefore, as the luminosity-temperature relation for groups which 
do not contain radio sources is unknown, we decided to investigate 
this by extracting radio-loud and radio-quiet subsamples from the 
group sample of Helsdon & Ponman. The criteria used in select- 
ing the radio-quiet subsample was that there were no radio sources 
associated with the group within 2 core radii of the central galaxy 
whose 1 .4-GHz luminosity was greater than 10 per cent of the value 
we obtained for NGC 3665, a radio galaxy (B2 1122-1-39) in the 
Helsdon & Ponman sample. We used NED, Simbad, FIRST and 
NVSS to check for radio sources. 

We compared this radio-quiet subsample with a sample of 
radio-loud groups consisting of the groups from the Helsdon & 
Ponman sample which were excluded from the radio-quiet subsam- 
ple, combined with the sample of Worrall & Birkinshaw (2000) 
(excluding 3C 449, as we have better constraints on its proper- 
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Figure 11. Plots of external pressure in the X-ray-emitting gas around 3C 66B (left) and 3C 449 (right) as a function of distance from the centre. Errors are 
la, determined from errors in the best-fitting /3-model parameters, which assumes an isothermal environment of 1.73 keV (3C 66B) and 1.15 keV (3C 449). 



ties from this work) which contains low-power radio-galaxy groups 
from the B2 sample. Fie. ll2l shows our two samples, with 3C 449 
and 3C 66B included for comparison. The best-fitting relation for 
each subsample is plotted. As there is a large intrinsic scatter in 
these relations, in order to determine realistic errors on parameters, 
we combined the measurement errors in temperature and luminos- 
ity with normally-distributed errors which are a small percentage 
of the temperature and luminosity. This extra error was increased 
until an acceptable fit was obtained for both samples; this required 
introducing an error of 0.005 x the logarithm of temperature and 
luminosity. The best-fitting relations are: 

log(Lx) = (4.42 ± 0.51) log(r) + (42.63 ± 0.09) (1) 

for the radio-quiet subsample (sample 1) and 

log(Lx) = (5.15 ± 0.52) log(r) + (42.27 ± 0.07) (2) 

for the combined radio-loud sample (sample 2), where Lx is the 
bolometric X-ray luminosity in erg s^^, and T is the temperature in 
keV. Both fits are consistent with a slope of ~ 4.9, the value found 
by Helsdon & Ponman; however, the intercepts for the two samples 
differ at greater than the 2a level. In addition we used a 1-D KS test 
to compare the temperature and luminosity distributions of sample 
1 and sample 2. We find that the luminosity distributions of the 
samples are consistent with their being drawn from a single parent 
population, whereas the temperature distributions are different at 
the ~ 90 per cent confidence level. 

Therefore, we find a small, but significant difference in the 
luminosity-temperature relation for groups containing radio galax- 
ies, with higher temperatures found at a given luminosity. There are 
two possible interpretations for this result: either the radio galaxies 
are heating their environments, or they are reducing the luminosity 
of the groups. To test whether the second scenario could be true, 
we assumed that all of the material through which the radio source 
passes while expanding has been affected by the radio source and 
no longer emits in X-rays (e.g. it may have been lifted out to a 
distance where the subsequent decrease in density produces this 
effect). Modelling the expanding radio source as a cone, and com- 
paring the square of number density integrated over the volume 
through which the radio source will have passed with the same in- 
tegral for the entire group, we find that only 7 per cent of the lu- 



minosity of 3C 66B could have been removed in this way. Using 
the original group luminosity, a temperature of 0.92±0.21 keV is 
predicted from our best-fitting relation for radio-quiet groups. The 
effect will be at a similar level in other sources. Therefore we con- 
clude that the radio galaxy could not have produced such a decrease 
in luminosity. 

For this reason, we interpret these results as providing indi- 
rect evidence for radio-source heating of group environments. This 
result shows that groups which contain detected radio sources are 
systematically hotter than those which do not, providing the first 
evidence that radio-source heating is occurring at a detectable level 
in the majority of such groups. It is apparent from Fig. ll2l that 3C 
66B is much hotter than predicted by either of the two relations 
plotted. There are three contributing factors which could explain 3C 
66B's exceptional behaviour Firstly, it is more powerful than any 
of the sources in the radio-loud sample, and so would be expected 
to transfer more energy to the group gas. As described above, it is 
also likely to be comparatively old for an FR-I radio galaxy, so will 
have been putting energy into its environment for a longer period 
of time. Finally, the environment of this source is less massive than 
the environments of other large, powerful sources, such as 3C 31. 
It seems plausible that these three factors have combined to make 
3C 66B particularly efficient at raising the temperature of its envi- 
ronment. 



5 CONCLUSIONS 

Our images of 3C 66B and 3C 449 provide direct evidence that 
their environments are having a major effect on the evolution of 
the radio sources. In both cases large rounded lobes are associated 
with large amounts of X-ray-emitting gas, whereas narrower lobes 
are seen where there is no surrounding material. In the case of 3C 
66B the presence of a blob of gas at the end of the eastern jet and 
bounded lobe is suggestive of an obstacle having a strong impact 
on the lobe morphology. 

We have revisited the luminosity-temperature relation for 
clusters and groups and find slightly different relations for sam- 
ples of radio-quiet and radio-loud groups. We interpret this as indi- 
rect evidence that radio-source heating is occurring in the majority 
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Figure 12. Plot of luminosity vs. temperature for two samples of X-ray 
groups. On the left is shown sample 1 (indicated by hollow squares): the 
radio-quiet subsample extracted from Helsdon & Ponman (2000), and sam- 
ple 2 (indicated by + symbols), the combined radio-loud subsample of Hels- 
don & Ponman (2000) and radio-galaxy sample of Worrall & Birkinshaw 
(2000). Best-fitting luminosity-temperature relations for each sample are 
plotted, as decribed in the text, illustrating that the sample containing radio 
sources show an increased temperature for a given luminosity. 3C 66B and 
3C 449 are included for comparison and marked with circles. 



of groups containing radio galaxies, and explains the higher-than- 
predicted temperatures we observe for the environments of 3C 66B 
and 3C 449. Better constraints on the temperatures of radio-quiet 
and radio-loud groups are needed to confirm this result. We find 
that either these two sources are expanding supersonically, or they 
are significantly older than expected from spectral age estimates. 
Supersonic expansion is unlikely due to the appearances of the 
sources, therefore we favour the second explanation, which implies 
that the sources are old enough and sufficiently powerful to have 
heated their environments and produced the higher temperatures 
we observe. 

The lobes of low-power radio galaxies are either out of 
equipartition or supported by a pressure contribution from parti- 
cles which we cannot observe directly. The lack of detected X-ray 
emission from the radio lobes rules out entrained thermal material 
at the temperature of the environment. We can also constrain the 
magnetic field strength in the lobes because of the lack of X-ray in- 
verse Compton emission, from which we conclude that if a depar- 
ture from equipartition explains the apparent pressure imbalance, 
the lobes must be magnetically dominated. 
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